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ABSTRACT 

Recent observations of the cosmic microwave background (CMB) have ex- 
tended measured the power spectrum to higher multipoles / >1000, and there 
appears to be possible evidence for excess power on small angular scales. The 
primordial magnetic field (PMF) can strongly affect the CMB power spectrum 
and the formation of large scale structure. In this paper, we calculate the CMB 
temperature anisotropics generated by including a power-law magnetic field at 
the photon last scattering surface (PLSS). We then deduce an upper limit on 
the primordial magnetic field based upon our theoretical analysis of the power 
excess on small angular scales. We have taken into account several important 
effects such as reionization and the modified matter sound speed in the pres- 
ence of a magnetic field. An upper limit to the field strength of \B\\ < 4.7 nG 
at the present scale of 1 Mpc is deduced. This is obtained by comparing the 
calculated theoretical result including the Sunyaev-Zeldovich (SZ) effect with 
recent observed data on the small scale CMB anisotropies from the Wilkinson 
Microwave Anisotropy Probe (WMAP), the Cosmic Background Imager (CBI) 
and the Arcminute Cosmology Bolometer Array Receiver (ACBAR). We discuss 
several possible mechanisms for the generation and evolution of the PMF. 

Subject headings: cosmology: cosmic microwave background — methods: numer- 
ical — magnetic field 
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1. Introduction 

The possible existence of a primordial magnetic field (PMF) is an important question in 
modern cosmology. The PMF could influence a variety of phenomena in the early universe 
(Tsagas and Maartens 2000; Grasso and Rubinstein 2001; Betschart et al. 2003; Bruni et al. 
2003; Clarkson et al 2003; Giovannini 2004; Vallee 2004; Challinor 2005; Kahniashvili et al. 
2005; Dolgov 2005; Kosowsky 2005; Gopal and Sethi 2005; Tashiro et al. 2005; Tashiro and 
Sugiyama 2005). However, there is still no firm understanding of the origin and evolution of 
the PMF, especially the observable magnetic field of 0.1 — 1.0 fi G in galaxy clusters (Clarke 
et al 2001; Xu et al. 2006). Recently, the origin of the PMF on such scales has been studied by 
a number of authors (Quashnock 1989; Boyanovsky et al. 2002; Bamba and Yokoyama 2004; 
Berezhiani and Dolgov 2004; Hanayama et al. 2005; Takahashi et al. 2005; Ichiki et al. 2006). 
There is, however, a possible discrepancy between theory and observation on the scale of 
galactic clusters. Temperature and polarization anisotropics in the CMB provide very precise 
information on the physical processes in operation during the early universe. However, the 
CMB data from the Wilkinson Microwave Anisotropy Probe (WMAP; Bennett, et al. 2003), 
the Arcminute Cosmology Bolometer Array Receiver (ACBAR; Kuo et al. 2004), and the 
Cosmic Background Imager (CBI; Mason et al. 2003) have indicated a potential discrepancy 
between the best-fit cosmological model to the WMAP data and observations at higher 
multipoles I > 1000. A straightforward extension of the fit to the WMAP data predicts a 
rapidly declining power spectrum in the large multipole range due to the finite thickness 
of the photon last scattering surface (PLSS). However, the ACBAR and CBI experiments 
indicate continued power up to / ~ 4000. This discrepancy is difficult to account for by a 
simple retuning of cosmological parameters. 

One possible interpretation of such excess power at high multipoles is a manifestation of 
the re-scattering of CMB photons by hot electrons in clusters known as the thermal SZ effect 
(Sunyaev and Zeldovich 1980). Although there can be no doubt that some contribution from 
the SZ effect exists in the observed angular power spectrum, it has not yet been established 
conclusively that this is the only possible interpretation (Aghanim et al. 2001) of the small 
scale power. Indeed, the best value of the matter fluctuation amplitude, as to fit the excess 
power at high multipoles is near the upper end of the range of the values deduced by the 
other independent methods (Bond et al. 2005; Komatsu and Seljak 2002). Toward a solution 
of this problem, it has recently been proposed (Mathews et al. 2004) that a bump feature 
in the primordial spectrum may provide a better explanation for both the CMB and matter 
power spectra at small scales. 

In this paper we consider another possibility. Inhomogeneous cosmological magnetic 
field generated before the CMB last-scattering epoch is also a plausible mechanism to pro- 
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duce excess power at high multipoles. Such a field excites an Alfven-wave mode in the 
baryon-photon plasma in the early universe and induces small rotational velocity pertur- 
bations. Since this mode can survive on scales below those at which Silk damping occurs 
during recombination (Jedamzik et al. 1998; Subramanian and Barrow 1998a), it could be 
a new source of the CMB anisotropies on small angular scales. Analytic expressions for the 
temperature and polarization angular power spectra on rather larger angular scales (/ < 500) 
by Mack et al. (2002) based upon the thin PLSS approximation have been derived for both 
vector and tensor modes. Subramanian and Barrow (2002) considered the vector perturba- 
tions in the opposite limit of smaller angular scales. Non-gaussianity in fluctuations from the 
PMF has also been considered (Brown and Crittenden 2005). A strong magnetic field in the 
early universe, however, may conflict with the cosmological observations currently available. 
The combination of those studies and current observations places a bound on the strength 
of the PMF of B < 1.0 - 10 nG. 

In order to compare the CMB anisotropy induced by the PMF with observations more 
precisely, we need to perform numerical calculations of the fully linearized MHD equations 
along with a realistic recombination history of the universe. In particular, we first need to 
develop a numerical method to predict the theoretical spectrum for intermediate angular 
scales, in which the analytic approximation becomes inappropriate. Recently, efforts along 
this line have been done (Lewis 2004) to obtain numerical estimates of the effects of the 
PMF on the CMB. Also Yamazaki et al. (2005a) have proven numerically that the PMF 
is an important physical process for higher multipoles I. They have constrained an upper 
limit to the PMF of B < 3.9 nG, based upon a likelihood analysis of the data (Yamazaki et 
al. 2005b). 

The WMAP data is very precise, but it does not extend to I >900. The PMF, however, 
affects the CMB power spectrum for higher multipoles / (Yamazaki et al 2005b; Lewis 2004). 
Thus, one purpose of this paper is to place a new limit on the strength of the primordial 
magnetic field by comparing our calculated result with the WMAP, CBI and ACBAR data, 
the latter two of which have measured CMB anisotropies for higher multipoles I than WMAP. 
We will also discuss the consistency between our CMB constraints and the magnetic field in 
galactic clusters. 



2. Primordial Stochastic Magnetic Field 

The Photon-baryon fluid can be assumed as single fluid since these particles interact 
each other very tightly before last scattering of photons. Since the trajectory of plasma 
particles is bent by Lorentz forces in a magnetic field, photons are indirectly influenced by 
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the magnetic field through Thomson scattering. We consider that the primordial magnetic 
field is generated from some moments during the radiation-dominated ear. In the case of the 
fait Friedmann-Robertson- Walker (FRW) background cosmology for the linear perturbation, 
since we treat the energy density of the primordial magnetic field as a first order perturbation 
and a stiff source for time evolution, the all back reactions from the fluid can be discarded. 

The electromagnetic tensor has the usual form 
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where Ei and Bi are the electric and magnetic fields. Here we use natural units c 
The energy momentum tensor for electromagnetism is 
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The Maxwell stress tensor, a %k is derived from the space-space components of the electro- 
magnetic energy momentum tensor, 
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The conductivity of the primordial plasma is very large and the electric field is negligible, 
i.e. E ~ 0. We call this condition "frozen-in" (Mack et al. 2002). This is a very good 
approximation during the epochs of interest here. In this way we obtain the following 
equations, 
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2.1. Two-Point Correlation Function 



We assume that the PMF, Bo, deviates from a homogeneous and isotropic distribution 
in a statistically random way. The power spectrum of fluctuations from a homogeneous 
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and isotopic distribution can then be taken as a power-law P(k) oc k nB ( iMack et al. 2002l ) 
where k is the wave number, and the spectral index n B can be either positive or negative. 
Evaluating the two-point correlation function of the electromagnetic stress-energy tensor, we 
can obtain a good approximation to the the vector isotropic spectrum for k < kc (Mack et 
al. 2002)0, 
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where the vector mode Lorentz force, L(k), is given by L(k) = fcnW(k), E A = |B A | is the 
magnetic strength for the comoving mean scale derived by smoothing over a Gaussian sphere 
for radius A, (where A = 1 Mpc in this paper). In this equation kc is the cutoff wave number 
in the magnetic power spectrum defined by 
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where h is the Hubble parameter in units of 100 km s -1 Mpc" 1 , (Mack et al. 2002; Jedamzik, 
et al. 1998; Subramanian and Barrow 1998a). We obtain the explicit time dependence of 
the primordial magnetic stress as Il^^k) = (k) / a 4 since Il^^k) are affected from the 
primordial magnetic field quadratically. The scalar mode of the magnetic field stress tensor 
n 1 - ) is obtained via a similar calculation (Koh and Lee 2002) 
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1 This has a form similar to that used in Mack et al. (2002) who treated two kinds of Fourier transforms 
of different normalizations. Here, we systematically adopt the same normalization for the Fourier transform. 
This removes the some of the uncertainty in the numerical calculations (Lewis 2004). 
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2.2. Ionization ratio 

In previous studies the ionization ratio x e in the early universe has been assumed to 
obey a step function in time, i.e. x e — 1 before the PLSS and x e = after the PLSS. 
However, more accurate theoretical results require the use of the exact ionization ratio. We 
have therefore carried out a numerical estimate of the ionization ratio by using the program 
"RECFAST" (Seager et al. 1999). In this paper we take exact account of the ionization 
ratio x e (r), and can now rewrite the Lorentz force as 

L(r,k) = fcx e (r)n«(k) . (10) 

3. Perturbation Evolution Equations 

To obtain the scalar and vector perturbation evolution equations we write the perturbed 
Einstein equation as 

(G^ + 5G^ V ) = 87rG(T Atl/ [ FL ] + 5T^ V [fl] ) , (11) 

where T^rpn is the energy- momentum tensor of a perfect fluid. The Einstein equation can 
then be separated into background and perturbation equations. The background Einstein 
tensor is given by 

G, v = - \^ V R , (12) 

while the perturbation tensor is 

5G^ = 5R^-U\5R . (13) 
The perturbed Einstein equation is then 

5R^-^JR = 8nG5T^ [FL] . (14) 
The perfect fluid energy-momentum tensor has the form 

^[FL] = P[FL]<A +(P[FL]+P[FL])CW, , (15) 

where the subscript [FL] denotes the fluid, while = dx^ j sj—ds 2 is the four-velocity of the 
fluid relative to an observer in the frame in which the Einstein equation ()14p is solved. We 
define the pressure p and energy density p for the ideal fluid which is instantaneously observed 
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by a comoving observer at rest relative to the fluid. We define the velocity fluid v % = dx l /dr, 
we can treat v l as the same order perturbation as5p[FL] = P[fl] — P[fl], ^P[fl] = P[fl] — P[fl] 5 
We obtain the energy-momentum tensor for linear order perturbations as follows, 
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Since we wish to consider the affect of the PMF, we combine the electromagnetic and perfect 
fluid energy momentum tensors, 
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We use the Boltzmann equations for the flat FRW Universe which Hu and White (1997) 
derives as follows. The moments of the temperature fluctuations G)["^ is 



e 



0«' 



(m) 



(21 



(2Z + 3) m 



1) z_1 



where 



(/ 2 -m 2 )(/ 2 



and the source term S[ m ^ is 



5*, 



( f c ef 
o 





/ 2 



f c ^ 0) + ^ 

r c ^ 1} + V 



(20) 



(21) 







f,p(2) - p 



(22) 



In equations (|20|) and (1221) we have introduced the optical depth r c between epoch r and the 
present epoch tq. 

pro 

r c (r')dr' , (23) 



Also in Eq. ( l22il we have introduced the anisotropic nature of Compton scattering as follow, 
The polarization evolutions follows, 
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and 
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3.1. Scalar Mode 

The matter fluid is assumed that it is representable as a ideal fluid and neglect the 
anisotropic pressure perturbations. Only adiabatic perturbations is considered and the en- 
tropy perturbations is negligible The linearized perturbation equations which are obtained 
(Ma & Bertschinger, 1995; Hu & White, 1997; Koh & Lee, 2000) from the Einstein equations 
up to first order are: 



+ H{j) + 20) 



3H 2 ip + 3H<p + k 2 <p 

2^-tf 2 W^(0-^) 
a J 3 



AnGa 2 5T^ 
A7iGa 2 {p + p)v {0 \ 

%Ga 2 5Tl 

127rGa 2 (p + p) 



(27) 
(28) 

(29) 

(30) 



where 
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while and a are defined as 



(33) 
(34) 



and 



(35) 



denotes the traceless component of T % y The conservation of energy- momentum is a conse- 
quence of the Bianchi identitiy; 
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Including the energy momentum tensor for the PMF and writing the magnetic energy density 
as 

P[em] = -y-, (37) 
Eq. ( 1361) leads to the following equations in fc-space: 

5 = -(l + w) (V 0) - 30) - 3H [j- p - w\ S , (38) 
1 + w op 1 + w 

where w = p/ p. Here, we can cancel a [em] m Eq. fl39|) in the linear perturbation limit. Thus, 
in the motion and continuity equations for the scalar mode, we can just add the energy 
density and pressure of the PMF to the general energy density and pressure respectively. 
From equations 038|) and (13"9"|) we obtain the same form for the evolution equations of photons 
and baryons as that deduced in previous works (Ma & Bertschinger 1995; Hu & White 1997). 

5, = -i;f-40 , (40) 

^ 0) = e^-a^+k^ + an^rivf-v^) , (41) 
S b = -vl 0) +3<P , (42) 

;,(0) _ «(0) . o, 2<lf , 4 Ay „ /,( ) .,(0) 



% = -~ v b 



+ c 2 s k\ + ^an e a T {v^ - vf ] ) + k 2 ip , (43) 

where n e is the free electron density, <jt is the Thomson scattering cross section, and er 7 in 
the second term on the right hand side of equation (141]) is the shear stress of the photon 
with the PMF. In the presence of the PMF, the magnetic pressure should be included in 
the acoustic term. Therefore, the sound speed in the second term on the right hand side of 
Eq.(|33l) should be 

c\ = c 2 bs + x e (r)U^/Ana 4 p b , (44) 

where x e (r) is ionization ratio (we have solved x e by using the RECFAST code (13 ll )) and c bs 
is the baryon sound speed without the PMF (see Adams et al 1996). 



3.2. Vector Mode 

The evolution of the vector potential V(r, k) under the influence of the stochastic PMF 
can be written (Hu and White 1997; Mack et al. 2002) as 
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where the dot denotes a conformal time derivative, pi and 7Tj are the pressure and the 
anisotropic stress of the photons (i = 7) and neutrinos (i = u). Since the vector anisotropic 
stress of the fluid is negligible generally, it is omitted. 

The magnetic field affects the photon-baryon fluid dynamics via a Lorentz force term 
in the baryon Euler equations. Following Hu and White (1997), the Euler equations for 
neutrino, photon and baryon velocities, v u , v 7 , and Vf, are written as 

€ ] - v = -k f^esl , (46) 

-V + r c (v? ~ « 6 W ) = -k ( ) , (47) 
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where R = (p-y+p^)/(pb+Pb) — (4/3)(p 7 /p;,) is the inertial density ratio between baryons and 
photons, while 6^ and 0^ are quadrupole moments of the neutrino and photon angular 
distributions, respectively. These quantities are proportional to the anisotropic stress tensors. 
Equations (l46l) - (l48l) denote the vector equations of motion for the cosmic fluid, which arise 
from the conservation of energy-momentum. 



4. Temperature Power Spectra 

We have constructed a Boltzmann code for the vector mode by explicitly expanding 
CMBFAST (Seljak & Zaldarriaga 1997). The temperature transfer function of the CMB 
today for the scalar and vector mode is expressed as an integral along the line of sight for 
/ > 2(Hu and White,1997; Mack et al., 2002) 
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where the ji(x) are the spherical Bessel functions with x = fc(ro — r), Z?i = ^//(/ + l)/2, 
-D2 = \/3-Di, and is the polarization term (Hu and White 1997). Using 0j (to, fc), we 
can compute angular power spectrum of temperature anisotropies: 



c < = 7 J dkk -wfiT ' (52) 

This power spectrum of temperature anisotropies is to be added to that induced by the 
scalar density perturbations and compared with observational data. 



5. Result and Discussions 



We discuss the effects of scalar and vector modes of the PMF and also the SZ effect on 
the CMB in this section. Note that the current data for higher I are known to be insensitive 
to the tensor mode which we ignored in the present calculations. Figures 1 and 2 show the 
CMB power spectra generated by the PMF. The effects of the PMF become progressively 
more important at higher multipoles and eventually dominates over the scalar fluctuations 
for / > 1000. As a result, the potential discrepancy between the observed and theoretical 
CMB temperature anisotropies at higher multipoles I is remarkably relaxed by the effects of 
the presence of a PMF (Yamazaki et al 2005). In the following subsections we discuss details 
on the behavior of the effects of the PMF on the CMB temperature anisotropies. 



5.1. Scalar and Vector modes 



The left hand side of Figure 1 shows the temperature anisotropies of the CMB power 
spectra which are generated by the PMF, and the right side hand of Figure 1 shows the 
primary temperature anisotropies calculated in the A Cold Dark Matter (LCDM) model with 
and without the effects of the PMF. The scalar part of the primary temperature anisotropies 
consists of several different sources of metric perturbations including photon and neutrino 
energy density perturbations as well as the perturbation due to the PMF. Therefore, the 
scalar mode generated by the PMF, as displayed by the green dash-dotted curves on the 
left hand side of Fig. 1, is evaluated by subtracting the two calculated CMB power spectra 
with and without the PMF in the LCDM model. Note that there is no such complication 
for the vector part of the primary temperature anisotropies because it arises uniquely from 
the vector mode of the PMF. 
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Both scalar and vector modes become larger with increasing field strength \B\\ of the 
PMF, but the variation of the vector mode is much stronger than the scalar mode (see 
the right hand side of Fig. 1). The scalar mode for the PMF gets complicated as is to 
be discussed later in this subsection. However the net effect is basically caused by two 
contributions: First, the increase of the sound speed for additional magnetic pressure to 
reduce baryonic matter fluctuation amplitude; and second, the changes of the metric and 
the total energy density which offset the first effect. Since the effect of the PMF for the 
vector mode is just an increase of the perturbation by the Lorentz force, the vector-to-scalar 
mode ratio increases as the strength of the PMF \B\\ increases. 

The PMF affects both the energy density and pressure in the scalar mode. The influence 
on the energy density is the same as that on the baryon density. The magnetic pressure, 
however, increases the sound speed of the fluid as discussed below Eqs. (H2l and ( 143|) . This 
produces interesting effects on the CMB (Adams et al 1996; Tsagas and Maartens 2000; 
Yamazaki et al 2006). In order to understand the effects of the change of the sound speed, 
c s , boosted by the PMF let us solve Eqs. (l42p and (143]) in the WKB approximation 



This solution clearly shows two different influences of the change of c s on the baryonic matter 
fluctuation: First, the increase of the sound speed due to the magnetic pressure results in 
a decrease in the amplitude of the evolution of 5b because of the prefactor of Eq. (l53l) . This 
means that in the presence of a magnetic field, the plasma pressure increases by the repulsion 
of lines of magnetic force. Thus, gravitational collapse of the matter is delayed. Second, in- 
creasing the sound velocity makes the frequency of the baryonic fluctuation 5b larger because 
of the exponent of Eq. fl53|) . This effect of the PMF on density field makes different kinds of 
amplitude boosts which depend on k (Yamazaki et al. 2006). The PMF shifts the acoustic 
oscillation for the baryon fluid either up or down depending on the combination of the sound 
speed c s and wavenumber k. Note, however, that since the amplitude of the PMF is so 
small, the difference between the shifted and non-shifted period of the acoustic oscillations 
essentially only depends on the wavenumber k. This effect is strongly constrained by the 
observed data as displayed on the right hand side of Figure 1. 

When one compares the scalar and vector modes, the scalar mode makes only a minor 
contribution to the CMB except at lower amplitudes. However, since the magnetic field 
pressure delays the gravitational collapse of matter, its effect on the matter power spectrum 
for large scale structure (LSS) is very interesting. Since in the present article we are only 
interested in the small scale structure associated with larger multipoles /, further details on 
this point will be discussed elsewhere (Yamazaki et al. 2006). 
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5.2. Dependence on B\ 

The CMB power spectrum generated by a cosmological magnetic field is generally char- 
acterized by a broad peak at / ~1000~3000 for a reasonable range of -2< ub <0 as shown 
in Fig. 2. We display five cases of the CMB power spectra for B\=l, 2, 4, 8, and 10 nG in 
Fig.l. As mentioned above, for lower / and B\, the scalar mode dominates (bottom three 
cases for B\ =1, 2 and 4 nG of the left hand side of Fig.l). For higher / and B\, the vector 
mode dominates. The amplitude is shown to increase with increasing magnetic field strength 
B\ for I > 1000, where the vector mode dominates. Some lengthy mathematical derivation 
of the power spectrum of fluctuations leads to C\ oc B\ for n B <-3/2, and C\ oc £?* 4 /( ns+5 ) 
for ub >-3/2 for higher multipoles I (Fig. 2). 



5.3. Dependence on % 

Fig. [2] shows the ub dependence of the CMB. For tib <-3/2, the second term in Eq.flTj) 
dominates. This fc-dependence causes a strong angular dependence for the vector mode 
perturbation. This results in changing the peak position for different multipoles I depending 
upon the power-law spectral index ng. On the other hand, peaks remain at nearly the same 
multipole I ~ 2000 independently of % for % >-3/2. This is because the vector mode 
perturbation from the magnetic field tends to be that of white noise when the first term in 
Eq.([7j) dominates. 



6. Constraint on the PMF in a MCMC Analysis 

To set limits on the PMF we have evaluated likelihood functions for fits to the WMAP, 
ACABAR, and CBI data over a wide range of the parameters, B\ and ub, for a stochastic 
PMF along with the usual six cosmological parameters, h,Qbh 2 ,{l m h 2 ,n s , A s , and r c . The 
quantity h is the Hubble constant in units of 100 km s _1 Mpc -1 . The quantities VL^h 2 and 
fl m h 2 are the usual present baryon and cold dark matter closure parameters. The quantities 
n s and A s are the spectral index and the amplitude of the primordial scalar fluctuation, 
respectively. The parameter r c is the optical depth. We adopt a flat cosmology. To explore 
this parameter space, we have employed the Markov chain technique (Lewis 2002). We also 
take account of the SZ effect in our analysis. For that, we adopt a fixed prior of ag = 0.9 
(Spergel et at. 2003; Komatsu and Seljak 2002). 



14 



6.1. Degeneracy 



In our likelihood analysis of the magnetic field parameters, we continued the MCMC 
algorithm until the cosmological parameters well converged to the values listed in Table 1. 
The inferred parameter values are consistent with those deduced in other analyses (Spergel 
et al. 2002; Mason et al., 2003; Kuo et al., 2004). We did not find obvious degeneracies of the 
magnetic field parameters with other cosmological parameters. We understand this for the 
following reasons. The primordial magnetic field is constrained (Jedamzik et al. 1998; Mack 
et al. 2002) by the cutoff scale for the damping of Alfven waves as indicated by Eq. (JHJ). 
Since this cutoff scale is small compared with the multipoles for the currently available 
data at I <C 3000, the manifestation of this effect is remarkably seen as a monotonically 
increasing contribution from both scalar and vector perturbations in the multipole region 
higher than I ~500 (Yamazaki et al. 2005a). This sensitivity of the CMB power spectrum 
to the primordial magnetic field differs completely from those to the other cosmological 
parameters, which helps resolve the degeneracies among them. There is, however, a strong 
degeneracy between the magnetic field strength \B\\ and the power spectral index ub because 
those parameters affect the amplitude simultaneously [see Eqs. ([7]) and (J9])]. 



In Fig. [3] we show results of our MCMC analysis using the WMAP , ACBAR and CBI 
data in the two parameter plane |B^| vs. ub- The lcr(68%) and 2<r(95.4%) C.L. excluded 
regions are bounded above by the thick curves as shown. We could not find a lower boundary 
of the allowed region at the lcr or 2a confidence level. Note, however, that we find a very 
shallow minimum with a reduced x 2 ~ 1.08. From Fig. [3] we obtain an upper limit to the 
strength of the PMF at the la(95.4%) C.L. of 



This upper limit is particularly robust as we have considered all effects on the CMB anisotropies, 
i.e. the effect of the ionization ratio, the SZ effect, and the effects from both the scalar and 
vector modes on the magnetic field, in the present estimate of |Ba| and ns- 



6.2. 



Limits on the PMF 



B x \ < 7.7 nG at IMpc 



(54) 



7. 



PMF Generation and Evolution 



In this section we discuss a multiple generation and evolution scenario of the cosmolog- 
ical primordial magnetic field that is motivated by the results of the present study. 
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To begin with we adopt the following three constraint conditions: 

1. A PMF strength \B X \ <7.7nG(lcr) at IMpc as deduced above by applying the MCMC 
method to the WMAP, ACBAR, and CBI data. 

2. The magnetic field strength in galaxy clusters is 0.1/iG < I-BcgI < ljU G (Clarke et 
al. 2001; Xu et al. 2006). Hence, if the isotropic collapse is the only process which 
amplifies the magnetic field strength, the lower limit to the PMF is ~ 1 — 10 nG for 
the PLSS. 

3. The gravity wave constraint on the PMF from Caprini and Durrer (20020 The bier- 
bang nucleosynthesis of the light elements depends on a balance between the particle 
production rates and the expansion rate of the universe. Since the energy density of the 
gravity waves pew contributes to the total energy density. However, p G w is constrained 
so that the expansion rate of the universe does not spoil the agreement between the 
theoretical and observed light element abundance constraints for deuterium, 3 He, 4 He, 
and 7 Li (Mag giore, 2000). 

Figures [3] and H] summarize the various constraints on the PMF by these three conditions. 
The region bounded by the upper red-solid curve is constrained by condition (i) as indicated, 
the green-dash-dotted horizontal line corresponds to the lower limit to the PMF from con- 
dition (ii), and the Black-solid, sky-blue-dotted, and pink-dashed lines, respectively, are the 
upper limit of the produced PMF from big-bang nucleosynthesis, the electroweak transition, 
and the inflation epoch. Figure. H] shows the allowed or excluded regions according to these 
multiple constraints depending upon when the PMF was generated: 

(A) . Nucleosynthesis: I+II+III region 

1.0 nG < \B X \ < 4.7nG 
-3 < n B < -2.40 

(B) . Electroweak transition: II+III region 

1.0 nG < \B X \ < 4.7nG 
-3 < n B < -2.43 



2 A in Caprini and Durer (2004) is for a scale of 0.1 Mpc, however our A is for IMpc. Thus, inclinations 
of lines in Figs. [3] and [4] are smaller than Caprini and Durer (2004). 
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(C). Inflation: III region 

1.0 nG < \B X \ < 4.7nG 
-3 < n B < -2.76 

Obviously, the upper limits on both \B\\ and n B become more stringent if the PMF 
is produced during an earlier epoch. Moreover, the limits we deduce are the strongest 
constraints on the PMF that have yet been determines. However, we caution that the 
evolution of the generation of the PMF during the LSS epoch is not well understood. Thus, 
if there are other effective physical processes for the generation and evolution of the PMF 
during the formation of LSS, our lower limit of the PMF parameters may decrease. In order 
to constrain the PMF parameters accurately, we should study the PMF not before but after 
the PLSS. 



8. Conclusion 

For the first time we have studied scalar mode effects of the PMF on the CMB. We have 
confirmed numerically without approximation that the excess power in the CMB at higher I 
can be explained by the existence of a PMF. For the first time a likelihood analysis utilizing 
the WMAP, ACBAR and CBI data with a MCMC method has been applied to constrain 
the upper limit on the strength of the PMF to be 

\B X \< 1-lnG . 

We have also considered three conditions on the generation and evolution of the cosmological 
PMF: 1) our result; 2) the lower limit of the PMF from the magnetic field of galaxy clusters; 
and 3) the constraint on the PMF from gravity waves. COmbining these, we find the 
following concordance region for the the PMF parameters; 

1 nG < \B X \ < 4.7 nG , -3.0 < n B < -2.4 . 

The PMF also affects the formation of large scale structure. For example, magnetic 
pressure delays the gravitational collapse. It is thus very important to constrain the PMF 
as precisely as possible. If we combine our study and future plans to observe the CMB 
anisotropies and polarizations for higher multipoles I, e.g. via the Planck Surveyor, we will 
be able to constrain the PMF more accurately, and explain the evolution and generation of 
the magnetic field on galaxy cluster scales along with the formation of the LSS. 
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A. Tight Coupling Approximation 

Since the Thomson opacity is larger before recombination, photons and baryons are 
tightly coupled. If the Thomson drag terms in Eq. pTI) . (143]) . (1471) . and (148]) are too large, 
it is difficult to solve these equations. Therefore, we here derive approximation equations 
in this limit. The tight coupling approximation for the scalar mode was introduced by Ma 
and Rertschinger (1997). Here, therefore, we only need to introduce the vector mode in this 
appendix. 

Using equations (T47j) and fT48]) . we obtain the following equations for the photons and 
baryons in the tight-coupling approximation. 




»i"-»? , = rl< 1, + y»SI • (A2) 



Writing Vj as vjp + (tty — v^), we can rewrite fjA2[) 

- v? = - ( + - + ^*6« ) , (A3) 



and using equations (I All) and (IA2I) . we get 

„W _ w (i) - 1 ( + _ jS 1 )) + ^i-B (1) 4- L(1) 



(A4) 
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Differentiating Eq. (IA4I) and ignoring terms of more than second order in f c , we get 



a 2R 



+ 



a (1 + R) 
1 

(1 + i?)r c 
0(t7 2 ) . 



a 

-a V » 



Substituting Eq. §EM into Eq. (jXTJ) yields 



= - 

7 i? 



b a b 5 72 a\ Pb + Pb ) 



(A5) 



(A6) 
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Table 1. Calculated LCDM model parameters which best fit the WMAP data with the 
effect of primordial magnetic field taken into account. 



Parameter 


Mean and 68% C.L. 


Errors 95% C.L. Errors 


Q b h 2 

QcDAlh 2 

n s 
A s 

h 


0231 + 00011 
u.uzoi„ 0010 

u.y (u_ 026 

u.»uo_ 12 9 
n 1 26+ - 034 

U.1ZU_ 039 

n 797+0.037 
u - ' z ' -0.037 


+0.0021 

-0.0020 

+0.024 

-0.018 

+0.054 

-0.046 

+0.236 

-0.194 

+0.11 

-0.10 

+0.074 

-0.070 


Note. 


Calculated corresponding cosmic expansion ages 


are t /Gjr - 


= 13.35^(68% C.L/ 


) and 13.35^;^(95% C.L.). 




Fig. 1. — Effects of the PMF strength \B\\ on the CMB. The five figures on the left hand 
side show the CMB temperature anisotropics from the PMF. The red bold curves are the 
temperature anisotropics from vector and scalar modes of the PMF. The green dashed- 
dotted curves are for the scalar mode, the skyblue doted curves are for the vector mode, 
and the purple thin curves show the SZ effect. The upper parts of the figures show the 
combined effects on the CMB. The lower parts show the effects of the scalar mode obtained 
by subtracting the LCDM model without the PMF from the LCDM model with the PMF. 
The five figures on the right hand side show the primary temperature anisotropics of the 
CMB with and without the PMF. The black short-dashed curves show the LCDM model 
without the PMF. The red bold curves are the LCDM model with the PMF. The purple 
thin curves are the SZ effect. The dots with lightgreeen, gray, and skyblue error bars are 
the WMAP, ACBAR and CBI data, respectivitly. On all figures, the power spectral index 
of the PMF is set to rag =-1.5. Also, the strengths of the PMF on both sides of this figure 
are set to be \B\\ = 8.0, 6.0, 4.0, 2.0, and 1.0 nG from top to bottom as labeled. 
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Fig. 2. — Effects power spectral index tib of the PMF on the CMB. The five figures on the 
left hand side show the temperature anisotropics of the CMB from the PMF, and the red 
bold curves are the temperature anisotropies from the vector and scalar modes of the PMF. 
The purple thin curves show the SZ effect adopted here. The five figures on right hand side 
show the primary temperature anisotropies of the CMB with and without the PMF. The 
black short-dashed curves are the LCDM model without the PMF. The red bold curves are 
the LCDM model with the PMF. The purple thin curves are the SZ effect. The dots with 
lightgreeen, gray, and skyblue error bars are WMAP, ACBAR and CBI data, respect ivitly. 
On all figures, the strength of the PMF is set to be \B\\ = 4.0 nG. Also the power spectral 
index of the PMF on both sides of this figure are set to = 0, -0.5, -1.0, -1.5, and -2.0 
from top to bottom as labeled. 
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Fig. 3.— Results of the MCMC method constrained by the WMAP, ACBAR, and CBI data. 
Excluded and allowed regions at the lcx(68%) C.L. and 2<r(95.4%) C.L. are shown in the 
two parameter plane of |B^| vs. ub, where |B^| is the primordial magnetic field strength 
and ris is the power-law spectral index. Solid and dotted curves are for Ax 2 =2.3 and 
6.17, respectively. The green dash-dotted horizontal line displays the lower limit of the field 
strength B\ = 1 nG deduced for the galaxy cluster scale at the PLSS (Clarke, et al 2001; 
Xu, et al 2006). The black-solid, skyblue-dotted, and pink-dashed lines are the upper limit 
of the produced PMF at the big-bang nucleosynthesis, the electroweak transiton, and the 
inflation epochs, respectively (Caprini and Durrer 2002). 
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Fig. 4. — Excluded and allowed regions at the lcr(68%) CL. and 2<t(95.4%) CL. in the two 
parameter plane |B^| vs. rig obtained by the MCMC method with the WMAP, ACBAR, and 
CBI data. |Ba| is the primordial magnetic field strength and ub is the power-law spectral 
index. The solid and dotted curves are for A% 2 =2.3 and 6.17, respectively. The black- 
solid, skyblue-dotted, and pink-dashed lines are the upper limit of the produced PMF at the 
big-bang nucleosynthesis, the electroweak transition, and the inflation epochs, respectively 
(Caprini and Durrer 2002). If the PMF is produced at the epochs of big-bang nucleosynthesis, 
the electroweak transition, or inflation, respectively, the region of I+II+III, II+III, or III are 
allowed by these constraints on the PMF for the galaxy cluster scale at the PLSS and the 
MCMC method with WMAP, ACBAR, and CBI data. 



